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The effect of in-plane electric field on Landau level spacing, spin splitting energy, average spin 
polarization and average spin current in the bulk as well as at the edges of a two-dimensional electron 
system with Rashba spin-orbit coupling are presented here. The spin splitting energy for a particular 
magnetic field is found to be reduced by the external in-plane electric field. Unlike the case of a two- 
dimensional electron system without Rashba spin-orbit interaction, here the Landau level spacing 
is electric field dependent. This dependency becomes stronger at the edges in comparison to the 
bulk. The average spin polarization vector rotates anti-clockwise with the increase of electric field. 
The average spin current also gets influenced significantly by the application of the in-plane electric 
field. 
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I. INTRODUCTION 

The proposal of the possible application of the spin 
degree of freedom in electronic devices has emerged as 
a new field called spintionicsi^— . The experimental real- 
ization of the spin-orbit interaction which lifts the spin 
degeneracy even in absence of magnetic field has boosted 
research interest in this field. The Rashba spin-orbit in- 
teraction (RSOI) is the particular interest of this field 
as it's strength can be enhanced significantly by ap- 
plying a suitable gate voltage^. The origin of the 
RSOI is due to the lack of structural inversion sym- 
metry in the semiconductor heterostructuresiS.. The 
RSOI provides a proficient way to control the electron's 
spin degree of freedom. The electric field perpendic- 
ular to a two-dimensional electron system (2DES) can 
enhance the Rashba spin-orbit coupling^. The RSOI 
has manifested itself through some remarkable new fea- 
tures in different transport properties^r^l. Moreover, 
spin Hall effect (SHE) is one of the most important 
consequences of the RSOI in a two-dimensional elec- 
tron system (2DES). In SHE, electrons with opposite 
spin polarizations get scattered in transverse direction 
of the applied electric field even in absence of magnetic 
field. There have been several theoreticaliSiilii^ as well 
as experimentaU^^ works on SHE. The successful detec- 
tion of spin accumulatior>iSiil at the edges has made the 
study of spin current interesting. 

The effect of an uniform in-plane electric field on 
the Landau levels in a two-dimensional system has 
been always an interesting quantum mechanical problem. 
The 2DES formed at the semiconductor heterostructure 
junction, which is also called conventional 2DES and 
graphene monolayer are two important examples of two- 
dimensional systems. The spacing between two consecu- 
tive Landau levels of a conventional 2DES is independent 
of electric field. On the other hand, the Landau level 
spacing in a graphene monolayer depends on the trans- 
verse electric field and it can be collapsed under suitable 
strength of the electric fieldi^"— . These two results are 



obtained by exact analytical calculation. The Landau 
levels of the Rashba spin-orbit coupled 2DES show two 
different energy branches with unequal level spacing—. 
The analytical derivation of Landau levels under in-plane 
electric field seems to be impossible. Therefore, we solve 
this problem by exact numerical calculation. 

The quantum Hall effect has been remained an at- 
tractive topic for physicists. The edge states and edge 
charge currents play very important role in exploring de- 
tails about quantum Hall effect^i^— . The spatially well 
separated left- and right-moving edge states give rise to 
non-dissipative transport^i^. In presence of the spin- 
orbit interaction, the transport of the spin-polarized edge 
channels may provide rich physics. Several theoretical 
analysis of the spin edge states in a Rashba system has 
been also reported^^i^. The crucial role played by the 
edge states in magnetization-^ and spin polarization^-? in 
2DES with RSOI has been studied extensively. 

In this work, we present the effect of in-plane electric 
field on Landau level spacing (LLS), spin splitting energy, 
average spin polarization and average spin current in the 
bulk as well as at the edges numerically. We find the Lan- 
dau level spacing and the spin splitting energy decreases 
with increase of the electric field. The average spin polar- 
ization rotates anti-clockwise as we increase electric field. 
The average spin current is modified due to the presence 
of the in-plane electric field. 

This paper is organized as follows. In Sec. II, we for- 
mulate the numerical method to solve the Hamiltonian 
of a single electron confined in a finite size 2DES with 
RSOI in presence of crossed electric and magnetic fields. 
We also discuss the spin-split Landau levels. Landau level 
spacing, average spin polarization and average spin cur- 
rent of an infinite 2DES with RSOI in presence of mag- 
netic field. In Section HI, we present numerical results 
and discussion of the effect of the transverse electric field 
on the spin splitting energy, LLS, average spin polariza- 
tion and average spin current. However, as the energy 
becomes different at the edges in comparison to the bulk 
for a finite system, the study of above mentioned prop- 



erties will also be covered at the edge with electric field. 
We present summary of our work in Section IV. 



II. FORMALISM FOR NUMERICAL 
CALCULATION 

We consider a spin-orbit coupled 2DES in a xy plane 
of area L^ x Ly which is subjected to an in-plane electric 
filed E = Ei and magnetic field B = Bz. The single- 
electron Hamiltonian with electronic charge q = —e \s 
given by 



H = 



(P + eA) 

2m* 



+ V{x) + eEx 



CTO 



+ -rl^xiVy + eAy) - cFyPx] + -^sBa^, (1) 

where m* is the effective mass of an electron, ctq is the 2 x 
2 identity matrix, ai {i = x, y, z) are the Pauli matrices, a 
is the Rashba spin-orbit coupling constant and V{x) = 
for \x\ < Lx/2 and otherwise infinity is the hard- wall 
confining potential along the x direction. Also, g is the 
effective Lande g-factor and jib = eft/(2me) is the Bohr 
magneton with me is the free electron mass. Here, we 
have chosen the Landau gauge A = Bxj. 

Before presenting numerical results, we shall review the 
bulk Landau levels of a 2DES with RSOI in absence of the 
electric field. The energy spectrum in the bulk can be ob- 
tained analytically by diagonalizing the Hamiltonian^^. 
The RSOI mixes the two spin components. For s = 0, 
there is only one level, same as the lowest Landau level 
without SOI, with energy E^ = Eq = {huc — <7Mb-S)/2 
with iOc ~ eB/m*. The corresponding wave function is 



* 



0,fc. 



,(r) = 



r^lkylJ 



(/)o(a; + xo) 



(2) 



where xq = kyl^ with magnetic length scale /q = 
\Jhl{eE). For a given quantum number s > 1 there 
are two spin-split branches of energy levels, denoted by 
-|- and — with energies 



Ef = shjJc±^E^ -h sEJuo, 



(3) 



where Eq, = 1wi*c? jf?. The corresponding wave func- 
tion for -|- branch is 



*^fc„(r) 



'''-y f Dscj^s^x + xo) 



yJLyAs 



5(2; + .To) ) ' 



and for — branch is 



*.,fc„(r) 



^""^ ( 0,_i(x + xo) 



yi^ \^ -Ds(t>s{x + XQ) 

\ + Dl with Ds 

\/ sEafvjJcl\BQ + ■\/£'q + sEahujc] and (psi^) 



(4) 



(5) 



where A, 



(l/-\/\/7r2''s!^o)e~^ /'^'■oHs{x/lo) is the normalized 
harmonic oscillator wave function. Here, s is the Landau 
level index. 

When both the hard-wall confining potential and the 
electric field are included, the Hamiltonian given by Eq. 
^ has to be solved numerically. To solve the Hamilto- 
nian numerically we choose the following wave function 
^i[x,y) — {\ / ^jLy)e^^yy 4>{x) with the function 0(x) ex- 
panded in the basis of the infinite potential well as 



cj>{x) 



i:-{Z(^+t)} 



(6) 



Now, using the time-independent Schrodinger equation, 
we get the following matrix equations for the spinors: 



[ 



ft27r2/2 



9 



^jLB<yzB 



2m*Ll 2 



0. (7) 



Here, the matrix elements are given by 



Fir, — 



2aL, 



Mo)2 



sm ■ 



X sm 



[n(e + ^)]de, (8) 



/•ir/2 

2cm / ' 


(if 


= / sm 


5 / f^ 


Lx J-ir/2 


I V 


. ^*^lLl 1 


»7r/2 


''" ^^ J 


-V2 



f)} 



d0. 



sm 



{io + l)}io + Oo 



X sm 



{"(^n)} 



d9 



(9) 



(10) 



and 



Ku 



2eEL., 



f_\in{l{e+^-)}0sin{n{6 



(11) 



where 9 = ttx/Lx, Oq — nxo/Lx- 

We solve these equations numerically in a truncated 
Hilbert space by considering 400 x 400 matrix Hamilto- 
nian and confirmed that first thirty eigenvalues in the 
bulk with E = exactly match with the results obtained 
from the analytical expression. Also, the probability den- 
sity of the corresponding low-energy states are in excel- 
lent agreement with the analytical results. We have nu- 
merically checked that in presence of the electric field 
E, the cyclotron center xq in the bulk is displaced by 
xe — eE/{m*uJe) from right to left side. For low- lying 
Landau levels, the electric field induced displacement xe 
is almost independent of s. When tq ~ the states 
are far away from the boundary and called bulk states. 
When Xq ss zLLx/2 the states are close to the boundary 
and called edge states. 



We define the spin splitting energy for a given Landau 
level s a,s A, = Et - Er . In the bulk with E = 0, 



As = 2-\/i?Q + sEahuc- It implies that the spin splitting 
energy A^ is increasing with increase of s, a and B. The 
LLS between two consecutive Landau levels is defined as 
AEf = Ef_^^ - Ef. The total of the LLS for spin-up 
and spin-down branches is AEg = A_B+ -f AEj = 2hu!c- 
This is independent of a and Landau level index s. 

The components of the average spin polarization are 
defined as 



p/-'.(t. 



lu'- 



t 



Kkjr) a^KkJr), (12) 



where i = x,y, z and A = ± represents the spin-up and 
spin-down components. The y component of the aver- 
age spin polarization is always zero. In absence of the 
electric field, the average spin polarization components 
can be evaluated analytically. In the bulk, the average 
spin polarization of the ground state (s = 0) is along the 
+z axis. In the excited states (s > 1), the average spin 
polarization in the bulk is P/"^- = X{h/2){l~\Ds\'^)/As. 
For a given s, the spin polarization of + branch is anti- 
parallel to that of the — branch. 

The spin current density for a given Landau level s 
is given by J^if'(r) = [^s.yY Jy''^s,y: where Jy' are 
the spin current operators with i = x,y, z. We use 
the conventional definition of the spin current opera- 
tor as Jy'' = hvyCTx/'i, and Jy" = h{vyaz + (JzVy)/^- 
The y component of the velocity operator is given by 
Vy = {x + Xo)uJc<7q + j-cTx- The average spin current is 

defined as (J,^-,) = Ea -^'jf' = Ea / <PrJ^;i'{^)- In ab- 
sence of the electric field, the average x-component spin 
current carried by a given Landau level can be easily ob- 
tained and it is given by {J ft) — a, independent of the 
Landau level index s. On the other hand, z component 
of average spin current is ( Jf ^) = 0. A point to be noted 
here that average spin current is linearly dependent on a 
and independent of magnetic field strength. On the other 
hand, the total spin current in absence of magnetic field 
is cubic in o*^. Following Ref32. , one can easily show 
that even in presence of electric field the spin current 



Js^y" is proportional to P^'^'' 



s,y 



and 



2m' 



This kind of relation does not hold between J^'f- 
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FIG. 1. Plots of the spin splitting energy in the bulk versus 
electric field for few s values. 
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FIG. 2. Plots of the spin splitting energy at the edges versus 
electric field for few s values. 



B = 0.5 T. To describe the left and right edge states, 
we have taken ky = zLOAALx/Iq- Therefore, the left and 
right edges are far away from each other in real space. 
We have also taken fc„ = to describe the bulk states. 



III. NUMERICAL RESULT AND ANALYSIS 



For the numerical calculation, we use the following pa- 
rameters: system size L^ — Ly — 1200 nm, electron's ef- 
fective mass m* = 0.068?7ie, Rashba spin-orbit coupling 
constant a = O.lao with ao = 10~^^ eV-m. Selective 
numerical parameters have been chosen to satisfy the fol- 
lowing conditions: cyclotron radius Re ^ lo < Lx/2 and 
the shift of the cyclotron orbit's center due to the max- 
imum applied electric field (cc_E)max + a;o < L^.- Here, 
(2;_E)max = 154 um for maximum E = 5 x 10^ V/m and 



A. Spin splitting and Landau level spacing 

In the previous section, we have seen that the spin 
splitting energy and the LLS in the bulk can be con- 
trolled by tuning the magnetic field. Here, we will see 
these quantities can also be controlled by the transverse 
electric field. The spin splitting energy and LLS are plot- 
ted in units of huc- In Fig. 1, we show how A^ in the 
bulk varies with E. Figure 1 shows that the spin split- 
ting energy is diminishing with the increase of the applied 
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FIG. 3. Plots of the Landau level spacing in the bulk versus 
electric field for few s values. 
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FIG. 4. Plots of the Landau level spacing versus electric field 
for few s values. Here, dashed and sohd lines represent left 
and right edges. 



transverse electric field. It implies that electric field ef- 
fectively reduces the effect of the RSOI in presence of the 
perpendicular magnetic field. In Fig. 2, we plot A^ at the 
two edges versus electric field. Comparing Fig. 1 and Fig. 
2, it is seen that the spin splitting energy at the edges is 
large compared to that of the bulk region when E = 0. 
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FIG. 5. Plots of the Landau level spacing in the bulk versus 
magnetic field for the fixed electric field. For better visualiza- 
tion, each curve is magnified by a suitable constant number 
as indicated in the figure. 



The spin splitting energy at the left edge is increasing 
with increase of electric field. On the other hand, A^. 
at the right edge is decreasing with increase of E. This 
is due to the fact that cyclotron orbit shifts from right 
to left with increasing electric field. The states around 
the left edge shifts towards the left boundary where as 
the states around right edge moves towards the bulk re- 
gion. Thats why the spin splitting energy of the left 
and right edge states increases and decreases with electric 
field, respectively. Moreover, the spin-splitting energy at 
the edges changes almost linearly with E whereas Ag in 
the bulk decreases non-linear ly with E. 

In Fig. 3, we show how LLS for spin- up and spin- 
down branches in the bulk varies with E. The LLS for 
upper branches decreases where as for lower branches 
it increases with the electric field. This is in complete 
contrast to the case of a 2DES without RSOI where LLS 
does not depend on the transverse electric field. But, 
it is similar to the LLS in graphene case. Moreover, for 
strong enough electric field, the LLS for spin-up and spin- 
down branches tends to saturate to huc which is the same 
as in the absence of RSOI. The reduction in the spin 
splitting energy and saturation in LLS in the bulk can 
be understood from the fact that when electric field is 
sufficiently high the effective Rashba coupling becomes 
very weak. It is interesting to note that AEg = AE^ + 
AE~ is always 2hu!c even in presence of electric field. We 
also plot LLS at the edges versus electric field in Fig. 4. 

In Fig. 5, we show the LLS as a function of magnetic 
field when the transverse electric field is fixed. It shows 
that both A+ and Aj increase with magnetic field. The 



slope of A+ is slightly higher than that of the A^ 
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FIG. 6. Diagrams of the average spin polarization vectors at 
various locations for different electric field strength. 



B. Average spin polarization 



FIG. 7. Plots of the average spin current for x and z com- 
ponents versus electric field in the bulk. The average spin 
current is in units of I'ilOc.Lx- 



When a suitable electric field is applied, the average 
spin polarization vectors for up and down states start 
to rotate anti-clockwise. The amount of spin rotation is 
solely determined by the electric field and location. The 
electric field effect on the spin polarization vector at the 
left edge is smaller compared to the bulk and right edge. 
This is because of the ky dependent asymmetric energy 
spectrum in presence of electric field. The anti-clockwise 
rotation of the average spin-polarization vector due to 
the electric field can be observed experimentally by Kerr 
rotation method. 



The average spin polarization vector versus electric 
field at different locations of the center of the cyclotron 
orbit is shown in Fig. 6. The increasing order of rows are 
accompanied with the increasing electric field whereas 
the columns are associated with the three different lo- 
cations of the cyclotron center i.e. left edge, bulk and 
right edge. The arrow line in the upper-half and lower- 
half circle stands for spin-up and spin-down branches, 
respectively. 

The top panel shows the average spin polarization vec- 
tor in absence of the electric field. In the bulk, the av- 
erage spin polarization vector is solely along the ±2; di- 
rection for spin-up and spin-down branches, respectively. 
At the edges the spin polarization vector lies in the xz 
plane. The magnitudes of the components of the spin po- 
larization vector depend on the location of the electron. 
It is interesting to note that the spin polarization vectors 
at the left and rights edges for -I- branch and — branch 
are not anti-parallel to each other. This is due to strong 
spin splitting at the edges as seen in Fig. 2. 



C. Spin current 

Figure 7 shows that in the bulk the average x com- 
ponent spin current carried by the Landau level s ~ 2 
is destroyed rapidly from the maximum value a due to 
the electric field, whereas the average z component spin 
current starts to increase from zero. After certain elec- 
tric field strength, it starts to decrease slowly with the 
increase of the electric field. We have checked that the 
upper and lower component of the z-component spin cur- 
rent individually follows the relation between spin polar- 
ization and spin current i.e; with increasing x-component 
of polarization the z-component spin current increases. 
But the total z component spin current does not follow 
this after a certain electric field strength as shown in the 
figure. 

In Fig. 8, we show how the average x and z compo- 
nents of the spin current carried by s = 2 Landau level 
change as we increase magnetic field for a fixed electric 



field. The spin polarization along z axis increases with 
increase of the magnetic field. In other words, the spin 
polarization along x axis decreases with increase of the 
magnetic field. The average spin current J^'^" will de- 
crease with magnetic field since Jf'if" = - g^A^Bg pA.s^ 
Our numerical result is consistent with the exact analyt- 
ical results. 




is under a perpendicular magnetic field. The spin split- 
ting energy in the bulk is diminishing with the applied 
transverse electric field. The Landau level spacing be- 
tween two successive Landau levels for upper or lower 
branches gets influenced by electric fleld. This result is 
in contrast to the 2DES without RSOI, but it is simi- 
lar to the graphene case. The LLS for upper or lower 
branches tends to saturate to hiOc for strong enough elec- 
tric fleld. These results indicate that strong transverse 
electric fleld can effectively reduces the effect of the RSOI 
in presence of magnetic fleld. The electric field also has a 
strong influence on the spin polarization vector depend- 
ing on the different location of electron. The average spin 
polarization vector rotates anti-clockwise as we increase 
electric field. The x-componcnt of the spin current is 
destroyed very fast by the applied electric field whereas 
the z-component of spin current increases initially then 
decays slowly. On the other hand, when electric fleld is 
fixed and magnetic fleld is varying, the 2;-component of 
the spin current decreases slowly after a certain magnetic 
field and x-componcnt of the spin current increases very 
fast. 



FIG. 8. Plots of the average spin current for x and z com- 
ponents versus magnetic field in the bulk. The average spin 
current is in units of eSf-Lx ■ 
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